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High-accuracy transmission measure- 
ments at an optical wavelength of 633 
nm and mechanical measurements of the 
thickness of a 13-fim thick silicon-crys- 
tal film have been used to calculate the 
absorption and extinction coefficients of 
silicon at 633 nm. The results are 
3105+62 cm-' and 0.01564±0.00031, 
respectively. These results are about 
15% less than current handbook data 
for the same quantities, but are in good 
agreement with a recent fit to one set of 
data described in the literature. 



Key words: absorption coefficient; etch 
stop; extinction coefficient; HeNe; high 
accuracy; silicon. 



Accepted: July 6, 1990 



1. Introduction 



The absorption coefficient [1] a at wavelength X 
is related to k, the imaginary part of the index of 
refraction, also called the extinction coefficient, by 



a=4 7r k/\. 



(1) 



where X is the vacuum wavelength of the incident 
radiation, which is related to the photon energy h v 
of the radiation by 



' Present address: Hughes Danbury Optical Systems, 
Wooster Heights Road, Danbury, CT 06810. 



100 



X A v= 1.23985 eV ju,m. 



(2) 



For silicon, k spectra have been derived from a 
Kramers-Kronig analysis of reflectance data [2], by 
inversion of ellipsometric data [3], and from a 
spectra derived from transmission measurements 
[4,5]. Neither ellipsometry nor a Kramers-Kronig 
analysis is well suited to the determination of small 
values oi k, so for silicon most k values quoted 
below about 2.5 eV have been determined from 
transmission measurements. Since the differences 
among the various data sets are fairly large, Taft [6] 
measured the absorption coefficient at the Hg 
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546.1-nm line (2.27 eV) for use in reducing ellipso- 
metric data on oxide films on silicon wafers. How- 
ever, at this time there is no consensus as to the 
correct spectra. This is illustrated by the fact that 
two recent compilations contain spectra that differ 
by much more than 20% at some wavelengths 
[7,8]. 

We have carried out a measurement similar to 
that of Taft at the 633-nm HeNe laser line. Our 
purpose was twofold. First, we wanted a smaller 
uncertainty at 633 nm than could be obtained from 
any of the published or compiled data. Second, we 
wanted to test the use of 1) an amplitude-stabilized 
laser as a radiation source, 2) extremely linear sili- 
con photodiodes as the radiation detector, 3) re- 
cently perfected etching techniques for sample 
preparation, and 4) state-of-the-art mechanical 
measurements of sample thickness. 

The advantage of using amplitude-stabilized 
lasers in conjunction with extremely linear photo- 
diodes is that it allows very precise transmission 
measurements to be carried out over a very large 
dynamic range with very high wavelength accu- 
racy and resolution. This was demonstrated re- 
cently in connection with a search for evidence of 
the second indirect gap in the silicon absorption 
coefficient spectrum [9]. 

The work reported here supplements the high- 
precision work reported in reference [9] by show- 
ing what additional steps are necessary to achieve 
high accuracy. Combining the techniques de- 
scribed here with those described in reference [9] 
should allow high-precision, high-accuracy mea- 
surements over a large range of transmittance. 

There are two reasons why a large range of 
transmittance is desirable. Since the absorption co- 
efficient of silicon changes by many orders of mag- 
nitude over the region of the indirect transition, the 
dynamic range over which accurate transmittances 
can be measured determines the number of samples 
needed to cover the whole spectral range. Second, 
there are certain advantages that accrue directly 
from calculating the absorption coefficient from 
small values of transmittance. The transmission of 
an optically thick sample is given by the simple 
expression 



/=(l-?-f) (l-?-b) exp(-ax). 



(3) 



where r^ and ?-b are the reflectances of the front and 
back surfaces of the sample, respectively, and x is 
the sample thickness. Therefore, 



da 
a 



6x 

X 



1 
ajc 



dt drf dry, 

t -^ {\~ri)-^ (\-r,)\ 



(4) 



Thus, the relative uncertainty in a will have its 
minimum value at the largest value of a x (smallest 
transmittance) for which the relative uncertainties 
in the transmittance and reflectances are indepen- 
dent of the transmittance. 

In the following sections, we describe the silicon 
sample, the optical measurements, the mechanical 
measurements, the data reduction, and the uncer- 
tainty analysis that were used to derive the absorp- 
tion and extinction coefficient values reported 
here. 



2. Silicon Sample 

The sample whose transmittance we measured 
was a 2-cm by 2-cm by 1 3-]am free standing film of 
silicon that was prepared as follows. A uniform 
thickness, p-type, epitaxial layer doped with about 
10'* B atoms per cm' (50 O-cm resistivity) was 
grown to a nominal thickness of 15 jam on a 100- 
mm diameter p""" (0.01 fl-cm) silicon wafer with 
(100) surface orientation. A nominal 25- by 25-mm 
square sample was cut from the wafer, and masked 
with wax so as to leave the central 2- by 2-cm area 
of heavily doped substrate exposed. The masked 
sample was then mounted on a sapphire carrier 
disk and etched in a rotating-drum mixture of 1 
part HF, 3 parts nitric acid, and 8 parts acetic acid 
[10]. After about 8 h, it was removed and cleaned. 
The result was a free-standing film, approximately 
15-ju,m thick, suspended in a 0.64-mm thick, 1.9-mm 
wide silicon frame. 

The sample was glued onto a flexible piece of 
plastic over a 2- by 2-cm hole cut in the plastic so 
that the etched side faced the plastic. The plastic 
served as a convenient handle for aligning and 
mounting the sample, and as a strain relief for the 
sample during handling. Although the sample was 
reasonably robust due to the frame, it still could be 
damaged by rough handling, and the plastic served 
to minimize stress on it during handling. 



3. Experimental Measurements 
3.1 Absolute Transmission and Reflection 
Measurements 

The absolute transmission measurements were 
carried out using the experimental setup shown in 
figure 1. A HeNe laser emitting 633-nm radiation 
was used as the source of radiation. The beam was 
passed through an electro-optical modulator and 
focused with a microscope objective onto a pin- 
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Figure 1. Apparatus for measuring the transmittance and reflec- 
tance of the silicon sample at the 633-nm HeNe laser line. 

hole aperture in an opaque screen. The radiation 
diffracted through the pin hole was collimated 
with a beam-expanding telescope to produce an 
Airy disk about 6 mm in diameter. An opaque 
screen with a 6-mm diameter hole (Airy Filter) was 
aligned with the first dark ring in the Airy pattern 
so that only the Airy disk passed through the aper- 
ture. This provided a fairly clean and uniform 
beam of about 40 /xW of 633-nm radiation. A beam 
splitter reflected some of this radiation to a photo- 
diode whose output controlled the transmittance of 
the modulator so as to maintain a constant power in 
the beam to within about 0.1% over a period of 
hours [II]. 

The amplitude-stabilized laser beam passed 
through the sample-holder section of a computer- 
controlled translation stage and was incident on a 
photodiode having a 1- by 1-cm active area that 
was located in position T in figure 1. The pho- 
tocurrent Iq from this photodiode was measured 
with a high-accuracy digital voltmeter and a trans- 
impedance amplifier. In this and all other pho- 
tocurrent measurements, the measured value was 
corrected by subtracting the dark current that was 
measured with the laser light shuttered and the rest 
of the apparatus in the same configuration as dur- 
ing the photocurrent measurement. 

3.1.1 Direct Transmittance and Reflectance 
The plastic sheet to which the sample was attached 
was positioned on the sample holder such that the 
laser radiation was incident on the etched side of 



the sample, and the bottom and top sides of the 
sample were parallel to one, and perpendicular to 
another of the directions of motion of the transla- 
tion stage on which the holder was mounted. The 
plastic sheet was fastened to the sample holder in 
this position with removable transparent tape. In 
this position the photodiode was 29 cm from the 
sample. The stages were then translated so that the 
6-mm diameter HeNe radiation beam was centered 
on the nominal center of the sample, and the pho- 
tocurrent Ii(J,k),j~ — 1,0, 1, A: = — 1, 0, 1, due to 
the directly transmitted radiation, was recorded as 
a function of the position of the sample over a 3 by 
3 grid of points on the sample having a 2-mm grid 
spacing, with the center of the grid pattern coinci- 
dent with the nominal center of the sample. 

The plastic sheet with the sample was removed 
from the sample holder and remounted as de- 
scribed above, but with the laser radiation incident 
on the polished surface, and the translation stage 
rotated about 9° off of the perpendicular to the di- 
rection of the laser beam to allow the reflectance of 
the sample to be measured by the photodiode in 
position R in figure 1. The resulting photocurrent 
/, reflected from the nominal center of the sample 
was recorded. The photodiode was relocated to 
position T in figure 1, but at 12 instead of 29 cm 
from the sample, and the photocurrent / from the 
radiation directly transmitted through the center of 
the sample was measured. 

The photocurrent / was 1.7% greater than 
7j(0,0), which was measured at the nominal center 
of the sample when the etched side of the sample 
faced the laser and the photodiode was 29 cm from 
the sample. This difference could be caused by a 
failure to irradiate exactly the same area of the sam- 
ple when aligning the laser beam with the nominal 
center of the sample, or by the different angle of 
incidence for the different measurements. It could 
also be caused by slightly more scattered radiation 
falling on the photodiode when located 12 cm from 
the sample than when located 29 cm from the sam- 
ple. In any event, the difference between / and 
7^(0,0) was not due to laser drift. After the mea- 
surement described above, the sample was re- 
moved from the sample holder and /o was 
remeasured. The two measurements of /o agreed to 
within 0.1%, indicating the stability of the ampli- 
tude-stabilized laser beam. 

3.1.2 Scattered Transmittance A scattermeter 
package [12] consisting of a box-shaped array of 
photodiodes with a 1- by 1-cm hole in the center of 
the array as shown in figure 2 was used to measure 
the transmitted radiation that was scattered out of 
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Figure 2. Illustration of the box-shaped scattermeter made from 
an array of photodiodes to measure the scattered component of 
the transmitted HeNe laser radiation. 

the directly transmitted beam. First, the photodi- 
ode used to measure the directly transmitted radia- 
tion was removed. Then, with the sample still out 
of the sample holder, the scattermeter was aligned 
in the beam so that it collected all of the beam on 
its interior surface, and the photocurrent /so was 
recorded. The sample was then remounted, as de- 
scribed above, with the etched surface facing the 
scattermeter. The scattermeter was then located 10 
cm from the sample as shown in figure 3, so that 
the 1- by 1-cm opening in the top of the scatterme- 
ter was in the position previously occupied by the 
photodiode measuring the directly transmitted ra- 
diation. The photocurrent /si was recorded. The 
scattermeter was then moved so that its base was 
only 2.2 cm from the sample, care being taken to 
keep the directly transmitted beam in the center of 
the 1- by 1-cm opening in the top of the scatterme- 
ter. The photocurrent /s2 was recorded. The re- 
sults, which are /si//o=0.060% and /s2//o= 0.094%, 
were used in the calculation of the total transmit- 
tance as described next. 

3.1.3 Transmittance and Reflectance Data Re- 
duction The transmittance t{j,k), j= — \, 0, 1, 
A: = — 1, 0, 1 of the sample at the same nine points 
on the 3 by 3 grid where the direct transmittance 
was measured was calculated as 



'"•*>=ll 



/^ , /sijr/s2 

h ho 



(5) 
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Figure 3. Illustration of the use of the scattermeter shown in 
figure 2 to measure some scattered components of the transmit- 
ted HeNe laser radiation. 

on the assumption that the transmittance of the 
sample was independent of whether the radiation 
was transmitted from the etched surface to the pol- 
ished surface, or vice versa, and that the transmit- 
ted radiation that was scattered scaled with the 
directly transmitted radiation. The t(j,k) values 
calculated from eq (5) are listed in table 1 . 

Table 1. The absolute transmittance t(J,k) measured on the sili- 
con sample described in this paper 



tU.k) 


»v 


-1 





1 


1 


0.00802 


0.00778 


0.00824 





0.00786 


0.00783 


0.00816 


-1 


0.00766 


0.00768 


0.00800 



The polished side of the sample, which faced the 
laser beam, was smooth enough to reflect specu- 
larly with negligible scattering. Moreover, when 
the radiation was incident on the polished side of 
the sample, the fraction that was reflected from the 
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etched side was attenuated by (l—rY exp(— 2 a 
x)<2X lO"* as it was transmitted through the pol- 
ished face of the sample to the etched face, and 
back again. Therefore, the scattered component of 
radiation reflected from the etched surface was 
negligible compared to the radiation reflected spec- 
ularly from the polished surface, and no attempt 
was made to measure it. 

The reflectances of both faces of the sample 
were calculated as 



rt=rf=///o=0.3454, 



(6) 



under the assumption that the roughness of the 
etched surface was not large enough to change the 
reflectance of the surface significantly, even 
though it was large enough to scatter a significant 
amount of radiation out of the directly transmitted 
beam. How this assumption was tested is described 
in the next section. 

3.2 Mechanical Measurements 

3.2,1 Surface Roughness The characterization 
of the roughness of the etched surface of the sam- 
ple was carried out by measuring the roughness 
along five lines located on the sample as shown in 
figure 4 using the NIST Stylus/Computerized Sys- 




Figure 4. Illustration of the location on the sample of the sur- 
face roughness measurements. Note that the length of the sur- 
face roughness scans is exaggerated in this figure. 



tem [13], in which a high-resolution TALYSTEP^ 
stylus instrument is calibrated on the 200,000 X 
vertical magnification range against a 0.029 l-jnm 
step-height standard that was measured interfero- 
metrically. To allow the fine structure of the sur- 
face to be resolved, a stylus width of 0. 15 ^um with 
a loading of 0.15 mgf (1.5 juN) was used. The 
TALYSTEP ISO filter (2CR) was used, and the 
cut-off wavelength was 0.0976 mm (25% attenua- 
tion). The traversing speed was 0.122 mm/s over a 
distance of 0.4 mm. Since the instrument's high-fre- 
quency response is 85 Hz, the surface profile could 
be detected over a wavelength range extending 
from about 1.44 jxm (corresponding to the fre- 
quency response of the instrument) to 97.6 fxm (the 
filter cut-off wavelength). The root-mean-square 
surface roughnesses R^ over this wavelength range 
varied from 7.2 to 10.0 nm for the five lines, while 
the arithmetic average roughness R^ varied from 
6.3 to 8.1 nm. This is not large enough to disturb 
the mechanical thickness measurements described 
next, nor to change the reflectance of the surface, 
but it is large enough to scatter both transmitted 
and reflected radiation significantly out of the re- 
spective direct and specular directions. 

3.2.2 Sample Thickness We used a high-accu- 
racy laser-interferometer micrometer (GCA 
LASERULER) having a resolution of 0.01 fxm and 
an uncertainty of ±0.13 jxm to measure the thick- 
ness of the sample. For these measurements, the 
micrometer was fitted with a 0.2-in (5.1-mm) di- 
ameter, hemispherical, steel-tipped probe, and a 
steel surface plate of nominal 130-mm diameter was 
centered under the probe. The digital output from 
the laser-interferometer that was connected to the 
probe shaft was set to zero when the probe was 
driven against the surface plate with the probe- 
driving mechanism. The probe was then raised 
from the surface plate and driven against it, and the 
zero-point height recorded six times. The probe 
was again raised, the sample placed on the surface 
plate with its polished surface down, and the probe 
driven to within a millimeter of the etched surface. 
In this position, the sample location was adjusted to 
be nominally centered under the probe as deter- 
mined by sighting along the diagonals of the sam- 
ple. 



^ Certain commercial equipment, instruments, or materials are 
identified in this paper to specify adequately the experimental 
procedure. Such identification does not imply recommendation 
or endorsement by the National Institute of Standards and 
Technology, nor does it imply that the materials or equipment 
identified are necessarily the best available for the purpose. 
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The probe was then lowered and raised six times, 
and the sample thickness was recorded each time. 
The zero-point thickness measurements were then 
repeated six times. The thickness h was calculated 
from the difference of the mean of the thickness 
measurements and the mean of zero-point measure- 
ments. The value determined in this way was 
h = 12. 86 ±0.02 }xm, where the stated uncertainty is 
one standard deviation of the reported value. The 
purpose of this value was to serve as a consistency 
check on the data on the variation of thickness 
over the sample that is described next. 

A small piece of graph paper was fastened to the 
surface plate with removable transparent tape in 
such a way that one corner of the plastic sheet to 
which the sample was fastened fell in the center of 
the graph paper when the sample was nominally 
centered under the micrometer probe. The graph 
paper did not extend under the sample for any posi- 
tion of the sample that allowed the probe to touch 
the etched surface of the sample. The plastic sheet 
was then aligned so that one edge was parallel to 
one of the lines on the graph paper, one corner was 
aligned with a point where horizontal and vertical 
lines on the graph paper crossed, and the microme- 
ter probe was located over the etched portion of 
the sample near one corner. The plastic sheet was 
then translated over an 8 by 8 grid with a 2-mm 
grid spacing and the height z(m,n), m =0, ..., 7, 
H =0, ..., 7, of the sample measured at each grid 
point. After the measurement of z(m,n) at each 
grid point, the plastic sheet was returned to the 
starting point (0,0) and the thickness z (0,0; m,n) 
was remeasured to monitor indirectly any zero 
drift in the digital readout of the micrometer. 

When the measurements described above were 
completed, the plastic sheet was removed and a 
zero-point reading was taken. The zero-point read- 
ing had drifted by 0.32 (xm during the set of 64 grid 
point and 64 reference-point measurements, but the 
reference-point measurements indicated that it had 
reached even larger values a few times during the 
measurements. Therefore, the reference-point mea- 
surements were used in making the zero-point cor- 
rections to the thickness map data. 

The plastic sheet was then alternately aligned at 
the point (1,3) and removed three times, so that 
three measurements of sample thickness followed 
by three measurements of the zero-point value 
could be recorded. The thickness ^(1,3) 
= 12.79±0.00 ^m of the sample at the point (1,3) 
was calculated as the mean of the sample thickness 
measurements minus the mean of the zero-point 
readings. The uncertainty was calculated in the 



same way as it was for the measurement at the 
nominal center of the sample that was described at 
the beginning of this section. 

Finally, the thickness h(m,n) of the sample at 
each grid point was calculated as 



h{m,n)=h(l,3)-]-[z{m,n)—z(0,0;m,n)] 
-[z(l,3)-z(0,0;l,3)]. 



(7) 



Table 2 shows the h{m,n) data obtained from eq 
(7). Figure 5 shows as crosses the approximate lo- 
cation of the 8 by 8 grid where the h(m,n) data 
were measured on the sample as determined by ob- 
serving the location of the sides of the micrometer 
probe relative to the corners of the sample after the 

Table 2. The thickness him, n) of the sample at points located at 
the approximate positions indicated in figure 5 

h{m,n) in ftm 



A' 



m 



1 



10.97 11.59 11.49 11.99 12.11 11.75 13.29 10.70 

12.42 12.66 12.80 12.58 12.60 12.56 12.24 11.68 

12.89 12.82 13.15 12.83 12.44 12.66 12.76 11.93 

12.94 13.01 12.79 12.81 12.54 12.80 12.88 11.84 
12.97 12.79 12.71 12.96 12.81 12.79 12.85 11.80 

12.95 13.05 12.95 12.66 13.02 12.88 12.50 11.68 



1 12.95 12.90 13.05 13.14 12.73 12.58 12.56 11.59 



12.17 12.61 12.68 12.74 12.55 12.30 11.97 11.36 
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Figure S. Approximate relative locations of the thickness mea- 
surements (+) and the transmittance measurements (x) relative 
to the center (circled X) of the sample. 
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measurement described above was completed. For 
comparison, figure 5 also shows as X 's the nominal 
location of the transmittance measurements re- 
ported earlier in this paper. The circled X is the 
point (0,0) for the transmittance data, and it is also 
the nominal location of the thickness measurement 
reported earlier in this section as a consistency 
check. That value, h = 12.86±0.02 /xm, agrees to 
within 0. 1 /xm with the expected value on the basis 
of the data shown in table 2. 

3,2.3 Elastic Compression Correction All of 
the silicon thicknesses reported so far contain an 
error due to the elastic compression of the hemi- 
spherical probe tip and the surface that it contacts 
under the load of a thickness or zero-reference 
measurement. The true thickness of the sample /;, is 
related to the difference h between the sample-in 
and sample-out measurement hy ht = h + Ah , where 
AA =a2— ai, as is illustrated in figure 6. To calcu- 
late Ah, we used the equation for the elastic com- 
pression a-, of a sphere in contact with a flat 
surface. 



( 3 77-/7 [Fflat+Fsphe 



.])' 



(8) 



where p is the total applied force, D is the diameter 
of the sphere, and 



F, = (l-0-/)/(77£,), 



(9) 



where ctj is Poisson's ratio of thejth material, and 
Ej is the modulus of elasticity of theyth material. 

Reference [14] gives 0.956x10"* inVlbf (1.39 
pmVN) for Fsteei- References [15] and [16] give 
0.44 and 1.56X10* Ibf/in^ (3.0 and 10.8 GN/m2) 
for CTsi and Esi, respectively, without indicating the 



crystallographic direction to which the values ac- 
tually apply. For forces applied to a (100) surface 
silicon, Esi is given by 

£si=(C„-F2 Cn) ■ (C„-Cn)/(C„ + Cn), (10) 

where Q; are the tensor stiffness constants of sili- 
con [17]. Reference [18] gives values for Q; from 
which we calculated £51=18.9X10* Ibf/in^ (130 
GN/m^). The diameter of the tip of the hemispher- 
ical probe was measured to be 0.20 in (5.1 mm), and 
the force applied during a thickness measurement 
was measured to be 0.154 Ibf (0.658 N). Using the 
above values, we calculated Ah =0.035 /xm and as- 
signed plus or minus the calculated value as a con- 
servative uncertainty. 



4. Data Reduction and Uncertainty 

Analysis 
4.1 Absorption Coefficient 

An absorption coefficient value was calculated 
for each of the transmittance measurements by in- 
verting eq (3) with rf=^b= 0.3454, as indicated in 
eq (6), to obtain 



aij,k) = 



where 



-ln{tU,k)/[(l-rd(l-r,)]} 
x(j,k) + Ah 



(11) 



x(j,k) = 

Zj X w(m—j — 3,n—k—3)h(m,n) 

m=j'+2 n=k-i-2 



108 



.(12) 




Figure 6. Illustration showing the relation between the true 
thickness A, of the sample and the measured thickness h due to 
the difference in compression between the sample-in and sam- 
ple-out measurements with the laser micrometer. 



fory = — 1, 0, 1, k = — l, 0, 1, and where the 
weights w(u,v) for « = — 1, 0, 1, 2, v= — 1, 0, 1, 2 
are listed in table 3. These weights were chosen as 
illustrated in figure 7 to average the measured sam- 
ple thickness over approximately the same 6-mm 
diameter area in the center of the sample as was 
used in the transmittance measurements. Table 4 
lists the resulting average thicknesses xij,k), and 
table 5 lists the absorption coefficient values ob- 
tained from eq (11). 

4.2 Uncertainty Analysis 

The uncertainty in the absorption coefficient can 
be calculated from the total differential in eq (4). 
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Table 3. The weights iv(u,v) used in calculating average thick- 
ness of the sample over the regions irradiated during the trans- 
mittance measurements 



V 



-1 



2 
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1 
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16 


14 


1 





11 


16 


16 
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-1 
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11 
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Table 4. The average thickness x(J,k) of the sample over the 
same areas for which the transmittance data are shown in table 1 



x<J,k) in nm 


*v 


-1 





1 


-1 


12.87 


12.75 


12.69 





12.84 


12.78 


12.79 


1 


12.89 


12.88 


12.81 



Table S. The absorption coefficients a(j,k) of silicon at 633 nm 
determined from the data in tables 1 and 4 

a(J,k) in cm"' 



.V 



-1 

3114 
3114 
3099 





3106 
3123 
3089 



1 

3090 
3128 
3081 



First, we consider the sources of error that are 
common to every point (/,^). These sources of er- 
ror include systematic radiometric errors in the di- 
rect and scattered transmittance measurements, and 
in the reflectance measurement, uncertainties asso- 
ciated with using the reflectance measured from 
the polished side of the sample to represent that 
from the etched side, and with using the sum of the 
measured scattered transmittances to represent the 
true scattered transmittance. Inaccuracies in the 
laser micrometer and the calculation of the depth 
Ah to which the micrometer probe penetrates into 
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Figure 7. Illustration of the averaging of the thickness measure- 
ments to approximate the same weighting over the 6-mm diame- 
ter region as obtained from the transmittance measurements. 
This weighting is appropriate for the center of the sample, but is 
a little skewed relative to other locations of the transmittance 
measurements on the sample. 

the silicon during the thickness measurement are 
also sources of error common to all of those mea- 
surements. These sources of error and our uncer- 
tainty estimates for them are listed in table 6. Also 
listed there is the quadrature sum of their contribu- 
tions to the relative uncertainty in the absorption 
coefficient. 

Referring to table 6 and eq (4), recall that the 
measured sample thickness was approximately 13 
ixm. Since the uncertainty associated with the laser 
micrometer is 0.13 /xm, it contributes 1% to the 
uncertainty in the absorption coefficient. Similarly, 

Table 6. The sources of error and associated uncertainty esti- 
mates for the sources of error that are common to the thickness 
and transmittance measurements at all points on the sample 



Source of error 



Estimated uncertainty 



dx/x 
Micrometer accuracy 
Compression correction 

de/e 

W,kyuo,o) 

I/Io 
/si/^so 

/S2/^S0 

(h^+h^yiso for true scatter 

d/T/(l-rf) 
V/o 

d/-b/(l-/^,) 
Use of rt for r^ 



1.00% 
0.27% 



1.0% 
1.0% 
1.0% 
1.0% 
3.8% 



<0.1% 



1.0% 



times 1 
1.00% 
0.27% 

times \/ax 
0.25% 
0.25% 
0.25% 
0.25% 
0.95% 

times 1/ou 
0.02% 

times 1/ou 
0.25% 



Sum in quadrature 



1.79% 
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the 0.035-ju,m uncertainty assigned to the compres- 
sion of the silicon, the surface plate, and the probe 
tip contributes 0.27%. 

There is a 1% uncertainty associated with the 
transmittance-photocurrent ratios because the gain 
of the current-to-voltage converter was 100 times 
greater for the measurements of the numerator 
photocurrents than for the denominator photocur- 
rents, and 1% is the between-range gain accuracy 
of the transimpedance amplifier used. An uncer- 
tainty of 25% of the measured components of the 
scattered transmittance was associated with the use 
of the sum of these quantities for the true scattered 
transmittance. Since the sum of the measured com- 
ponents of the scattered transmittance was about 
15% of the total measured transmittance, the rela- 
tive uncertainty is about 3.8%. An uncertainty of 
less than 0.1% is associated with the reflectance 
ratio since the same gain was used for both the 
numerator and the denominator photocurrents in 
this measurement. However, an uncertainty of 1% 
was associated with r^, to allow for the possibility 
that it was smaller than rj due to the roughness of 
the etched surface. According to eq (4), each of 
these uncertainties is multiplied by a factor of \/cxx 
(about 1/4 for this sample and wavelength) in its 
contribution to the uncertainty in the absorption 
coefficient, as shown in table 6. 

Sources of error that contribute differently to 
the absorption coefficient values at the different 
points (j,k) are the random errors in the measure- 
ments of the transmittance, reflectance, and thick- 
ness. The first two are negligible at the 0. 1 % level. 
The third might not be, but its actual contribution 
to our data cannot be distinguished from the error 
associated with the variation in thickness of the 
sample as a function of position over its surface. 

It is not really clear how to estimate the uncer- 
tainty to be associated with the variations in sample 
thickness. On the one hand, the thickness averag- 
ing procedure was designed to eliminate errors 
from this source. On the other hand, the procedure 
is not perfect because the transmittance values and 
the average thicknesses do not weight the same 
parts of the sample in the same way. The transmit- 
tance values weight the points of the sample within 
a 6-mm diameter region according to the intensity 
of the Airy disk. The average thicknesses weight 
13 points on a 2- by 2-mm grid nonuniformly to 
approximate equal weighting of all points in a 6- 
mm diameter region. Furthermore, the centers of 
the corresponding regions over which the transmit- 
tance and thickness averages were computed are 
only approximately aligned. 



4.3 Reported Value and Uncertainty 

After trying a number of different approaches to 
the uncertainty analysis, we decided to use the av- 
erage of the a(J,k) values in table 5 for the re- 
ported value, and to use the half-range of the 
values as the uncertainty. This results in a contribu- 
tion of 0.76% to the relative uncertainty in the ab- 
sorption coefficient. When we add this relative 
uncertainty in quadrature with the relative uncer- 
tainties reported in table 6, we obtain 1.94%. 
Therefore, we report a = 3105 cm"' and 
A =0.01564 with an estimated uncertainty of 2% 
for silicon at a vacuum wavelength of 633.00 nm or 
a photon energy of 1.9587 eV. 

4.4 Comparison with Existing Data 

Figure 8 compares our data at 1.96 eV with some 
of the more noteworthy data previously reported 
in the 1.9- to 2.3-eV spectral region. The open cir- 
cles are the handbook data of reference [7] and the 
closed circles are the handbook data of reference 
[8]. The dashed line is the data (for which no tabu- 
lated values were published) of reference [4], and 
the full line is data (for which no tabulated values 
were published) of reference [5]. The open dia- 
mond is the value reported in reference [6], along 





U.UiJ 




O 


1- 


0.04 


- 


o 


LU 








O 

LL 


0.03 


_ 


o o3" 

GO O >'^J^ 


LU 






o o y^ .^^ 


o 






rs •^ ^ 


o 






O ^ ^ 


2 






^-^ ^^ 


O 


0.02 


_ 


•■^ .y^ 


1— 




•o 


■^ ^ 


o 






X 








LU 


n ni 


i 


1 1 1 1 1 1 



1.9 2.0 2.1 2.2 2.3 

PHOTON ENERGY (eV) 

Figure 8. Comparison of Handbook data (open circles from ref- 
erence [7] and closed circles from reference [8]), published data 
reported without tabulation (dashed line from reference [4] and 
continuous line from reference [5]), and single-point transmit- 
tance measurement (open diamond from reference [6]) with the 
value reported at 633 nm in this paper. 
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with the reported uncertainty. Our value is the 
closed diamond. 

The result reported here has an estimated uncer- 
tainty that is significantly smaller than the other 
results shown in figure 8, due to combining a 
higher quality sample with more accurate and pre- 
cise thickness and transmittance measurements. 
Two sources of uncertainty dominate the uncer- 
tainties listed in table 6. These are the measurement 
of the thickness of the sample and the measurement 
of the scattered component of the transmittance. It 
should be possible to improve both of these mea- 
surements to reduce the overall uncertainty to 
below 1%, either by using more accurate instru- 
mentation, and/or by improving the sample qual- 
ity. In either case it would be necessary to obtain a 
more accurate laser micometer. Other required im- 
provements in instrumentation would include a 
more accurate scattermeter, and a precision trans- 
lator that clamps onto both the surface plate used 
in the thickness measurements and the sample 
holder used in the transmittance measurements. 
The latter would allow precise alignment of the 
areas weighted by the transmittance and thickness 
measurements, as well as a finer density of points 
for the thickness measurements. Improvements in 
the sample preparation procedure that resulted in a 
sample having a more uniform thickness and a 
smoother rear surface would eliminate the require- 
ment for the improved scattermeter and the trans- 
lation stage. However, at the current time, it is not 
clear how to obtain the necessary improvements in 
sample preparation. 
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